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We explore the model-independent constraints from cosmology on a dark-matter particle with 
no prominent standard model interactions that interacts and thermalizes with other particles in 
a hidden sector. Without specifying detailed hidden-sector particle physics, we characterize the 
relevant physics by the annihilation cross section, mass, and temperature ratio of the hidden to 
visible sectors. While encompassing the standard cold WIMP scenario, we do not require the freeze- 
out process to be nonrelativistic. Rather, freeze-out may also occur when dark matter particles 
are semirelativistic or relativistic. We solve the Boltzmann equation to find the conditions that 
hidden-sector dark matter accounts for the observed dark-matter density, satisfies the Tremaine- 
Gunn bound on dark-matter phase space density, and has a free-streaming length consistent with 
cosmological constraints on the matter power spectrum. We show that for masses < 1.5 keV no 
region of parameter space satisfies all these constraints. This is a gravitationally-mediated lower 
bound on the dark-matter mass for any model in which the primary component of daqrk matter 
once had efficient interactions — even if it has never been in equilibrium with the standard model. 



I. INTRODUCTION 

For many decades we have been aware that dark mat- 
ter exists through its gravitational effects on galaxies and 
the Universe as a whole. ^ Over this time, compeUing can- 
didates for dark matter have been proposed that often 
try to link the origin of dark matter to other outstand- 
ing problems in physics and cosmology. Amongst others 
these include, for example, weakly interacting massive 
particles (WIMPs) (e.g., [2]), axions [3], sterile neutri- 
nos [1], Super WIMPs [S], and asymmetric dark matter 
(ADM) models that link the origin of dark-matter and 
baryon number (e.g., [6j). Many of these models are al- 
ready well constrained from a remarkable effort in exper- 
imental dark matter physics including direct and indirect 
searches for dark matter and constraints from cosmology, 
and improved measurements and data from the Large 
Hadron Collidor (LHC) have the potential to further shed 
light on dark matter physics (see, e.g., Ref. [7])- 

Despite these intriguing possibilities for dark matter 
physics it remains entirely possible that dark matter be- 
longs to hidden sector and lacks standard model gauge 
charge. In such a case dark matter is expected to have at 
most extra-weak interactions with the standard model, 
and in the extreme case interacts only via gravitational 
forces. In such a case is dark matter physics wholly in- 
accessible to us? We take heart in the knowledge that 
the only reason wc know about dark matter is via its 
gravitational effects in cosmology and astrophysics. 

The most scrutinized candidate for dark matter is ar- 
guably the WIMP because the so-called "WIMP miracle" 
ensures that the simple process of its weak-scale interac- 
tions freezing out from a thermal plasma give it nearly 
the right relic abundance to be all of the dark matter. 
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But as this result only relies on the size of the total an- 
nihilation cross section being close to typical electroweak 
value (Tew ~ 10~® GeV^ (and not that the annihilation 
products are standard model particles) the WIMP mira- 
cle can be extended to a "WIMPless miracle" in hidden- 
sector models that also have a similar value for the cross 
section but with a range of masses extending down to 
0(10) keV [8 . In Ref. [5] a gauge mediated model of 
supersymmetry breaking is adopted that naturally gives 
a ~ (Tew and WIMP-like dark matter that decouples 
from the hidden plasma when it is nonrelativistic. Also 
see Ref. for different possible origins of hidden sector 
dark matter. In Ref. [10] hidden-sector models that de- 
coupled when ultrarelativistic were considered. In this 
work we present a unified treatment of the freeze-out of 
thermal relics in hidden sectors for arbitrary cross sec- 
tion (T, only requiring that constraints from cosmology 
are satisfied, and find viable dark matter that freezes out 
when it is relativistic, semirelativistic, or nonrelativistic. 

Already, there are stringent bounds on hidden sector 
dark matter from cosmology if the hidden-sector reheated 
to a temperature at or above the temperature of the vis- 
ible sector. If the hidden sector has more than a few 
light species for dark matter to annihilate into then con- 
straints on the number of relativistic degrees of freedom 
from big bang nucleosynthesis (BBN) rule out such mod- 
els [TTVH3j . However, in hidden sector models the reheat 
temperature of the hidden sector can easily be different 
from that of the visible (standard model) sector depend- 
ing on the details of the reheating process or if the visi- 
ble sector and hidden sector contain a different particle 
spectrum (different number of degrees of freedom as a 
function of temperature) and thus cool differently. Note 
that implicit in saying the visible and hidden sectors are 
at different temperatures is the assumption that the two 
sectors are not in thermal contact (or at least they lost 
thermal contact considerably before BBN) and processes 
that could thermalize them in a Hubble time are ineffi- 
cient. 

As discussed above, we take a model-independent ap- 
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proach to hidden-sector dark matter in this paper and 
let the cross section a be a free parameter. In addi- 
tion wc aUow the mass of our dark matter particle x 
and the hidden-to- visible temperature ratio ^ = T^/Tf 
(at dark matter freeze-out) be additional parameters and 
map of the viable region of this three-dimensional param- 
eter space. Depending on the values of m^, ^, and cr the 
freeze-out process can be relativistic, semirelativistic or 
nonrelativistic. 

Our paper is organized as follows: In Sec. [IT] we dis- 
cuss the standard Boltzmann equations involved in the 
freeze-out process and how they must be extended for 
the general hidden-sector case with ^ < 1. In Sec. |III| we 
discuss our approach to modelling the temperature de- 
pendence of the annihilation cross section in detail, how 
we treat the problematic semirelativistic regime, and a 
numerical approach to relativistic and nonrelativistic de- 



coupling within the same framework. In Sec. IV we calcu- 
late bounds arising from cosmological limits on the free- 
streaming length of dark matter, while in Sec.|V]we cal- 
culate limits from the generalized Tremaine-Gunn bound 
on the phase space density of dark matter. We discuss 
our results in Sec . |VH the strength of hidden-visible in- 
teractions in Sec. |VII[ and conclude in Sec. |VIII| 



II. 



THE RELIC ABUNDANCE 



In this section we review the standard Boltzmann 
equation for the freeze-out of annihilation of dark matter 
and its extension to the hidden-sector case. Suppose we 
have dark-matter particle of mass interacting with a 
plasma of other particles at temperature T. In terms of 



the yield, or number density to entropy ratio, Y 
we can write 



dY 1 s{m-^) 



M [Y^-YoixY' 



(1) 



where 



; m^/T a convenient time parameter, Yq{x) 
is the equilibrium yield (the yield if dark-matter annihi- 
lation were in equilibrium), and s{'m^) and H(m^) are 
the total entropy and Hubble rate evaluated at T = 
respectively.^ In the Maxwell-Boltzmann (MB) approxi- 
mation wc can write 1141 



Yo{x) 



45 d^ 
47r4 g^{x) 



x^K2ix) 



(2) 



where d-^ is the number of internal freedoms of a dark- 
matter particle, g^{x) is the effective number of en- 
tropy degrees of freedom, and here and elsewhere, Kn{x) 
are modified Bessel functions. This can be expressed 
in a convenient form by using the scaled variable y = 



[s{m^) / H (m-^)]{av)Y . In the limit that the temperature 
dependence of (av) is negligible we then have 



dy 
dx 



(3) 



where yo is the equivalently-scaled equilibrium yield. If 
dark matter is part of the visible sector and decouples 
when nonrelativistic then standard freeze-out results are 
found by solving Eq. Q with yo{x) as the initial condi- 
tion for y at small x (high temperatures). 

However, if dark matter is part of a hidden sector de- 
coupled from the visible sector when it freezes out then 
its equilibrium distribution is given by the hidden-sector 
temperature Th — ST rather than T, where we have in- 
troduced the hidden-to- visible temperature ratio ^. It's 
conventional to cast the total density (that appears in the 
Hubble rate) or entropy density in terms of an effective 
number of degrees of freedom g{T) [or entropy degrees 
of freedom g^iT)] at a temperature T. As discussed in, 
for example, Refs.pl fTO]. in ^ universe with both a visible 
and hidden sector there are two copies of these functions, 
gv{T) for the visible sector and gh{Th) for the hidden 
sector. Since for relativistic particles the energy density 
scales as p oc T"* and entropy density scales as s oc T'^ we 
find that the combined g{T) and g^{T) take the form 



5(T) = 5„(T) 
g^T)^gl{T) 



9v{T) 



(4) 



where the rightmost approximation typically holds when 
f < 0.5 and is used it in the remainder of the paper. If the 
hidden sector is relatively cold compared to the visible 
sector then the energy density (and Hubble expansion 
rate) of the Universe is essentially given by the density 
of visible-sector radiation. Thus the key change to Eq. ([T]) 
is just that Yq{x) — > yo(x;^), where 



Ya{x-i) = S 



45 



47r4 g''{x) 



x^K2{x/0 



(5) 



is the equilibrium yield of dark matter in a hidden sector 
and the standard case is recovered when ^ — >■ 1. In terms 
of the scaled yield y we have 



dy 
dx 



1 



y 



(6) 



with 2/0(2^; C) = [si'm^)/H{'m^)]{av)YQ{x;£,). Equation 
^ is valid when the thermally averaged cross section is 
almost constant at freeze-out. We discuss in the next 
section the more general case when the temperature de- 
pendence of the cross section cannot be neglected. 



III. THERMALLY AVERAGED CROSS 
SECTION 

While approximate analytical solution are well known 
when dark matter decouples when either highly relativis- 



^ We use units with c = 1 and ki 



tic Xf 



(/Tf ^ 1 or nonrelativistic x/ > 3, for the 
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intermediate semirelativistic case the thermally-averaged 
cross section (tru) cannot be expanded in either dark 
matter mass or velocity, and analytical methods are less 
tractable. However, in Ref. [TS], a useful and accurate 
interpolating expression for the thermal average cross 
section was introduced that can be used to track the 
dark matter freeze-out process even in the semirelativis- 
tic regime x/ ~ 1. The general expression for {av) is 
given by 



1 



{av) 



dscris-imDy/sKiiy/s/T) 



(7) 

While, strictly speaking, using the MB distribution is 
improper for the semirelativistic case, it has been shown 
to yield accurate results for the final relic density [IB] . 
This is partly due to cancellations between the numer- 
ator and denominator and also because the final relic 
density becomes quite insensitive to the precise value of 
Xf as freeze-out moves from the nonrelativistic toward 
the relativistic case. In the case of an s-wave process,"^ 
after some simplification, (av) can be put in the form 
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FIG. 1: The (scaled) yield y{x) and equilibrium yield y{x;S^), 
where x^is — x = m^/T represents visible sector tempera- 
ture, are shown for different choices of a and ^. The solid, 
dotted and dashed curves represent ^ = 0.9, ^ = 0.3, and 
^ — 0.05 respectively for a = O.lcro and ctq = 10~^ GeV~^. 
In agreement with Ref. [S], increasing ^ at fixed a increases 
the relic density. 



dt {f + x'^yKi{2x/t^ + x^) 



(8) 

Now in the nonrelativistic and ultrarclativistic limits it 
reads 



a;S>3 

M«^limM = ^ 



and 



;{12 + 5X^ 



(9) 



Xf = m^/Tf. Exactly how Tf is defined is somewhat 
arbitrary, but a natural definition is the temperature Tf 
at which T < H for the first time — the interaction rate 
drops below the Hubble rate. Such a definition can be 
used to accurately estimate the final yield in the relativis- 
tic and semirelativistic cases, and can naturally extend 
to the conventional nonrelativistic case. Explicitly, we 
use the condition 



repectively. A simple ansatz that interpolates between 
the these cases is 



[sY{av)] 



HI 



(12) 



(av) = (Jj{x) 



(10) 



where f{x) is chosen to smoothly reproduce both the ul- 
trarelatvistic (x <C 1) and nonrelativistic {x 3) cases. 
Incorporating this new temperature-dependent cross sec- 
tion {av) — af{x) wc find Eq. ^ takes the form 



dx 



[y^ -yo(2;;0] +v 



din fix) 
dx 



(11) 



With the above approximation, we can numerically solve 



Eq. (11) for the freeze-out process as a function of a, 
m^, and ^. A key quantity we need to know to es- 
timate the final relic abundance and whether it occurs 
in the ultrarclativistic, semirelativistic, or nonrelativistic 
regimes is the freeze-out temperature Tf or equivalcntly 



which in terms of the scaled yield becomes simply 
y{xf) = Xf. We use this later condition to determine Xf 
numerically. The freeze-out solution is shown in Fig. [T] 
for 3 different choices of ^ at fixed a, where we see higher 
f results in higher relic density. Despite higher values 
oi Xf, and thus later freeze-out with more exponential 
suppression, the larger relative temperature more than 
compensates in this example. 

Although reasonable approximations can be found in 
the nonrelativistic and ultrarclativistic cases, to calculate 
the present relic density in what follows we use the yield 
found from numerically solving Eq. (Ill at a late time x^o 
where the yield is essentially constant. The cosmological 
dark-matter density is then 



r2^ — m^SQ 



(13) 



^ The p-wave case can be treated similarly althougli tlic approxi- 
mating ansatz is different |15l . 



where pc is the critical density of the Universe. We use 
the measured value of the dark matter relic density from 
Ref. |16) to constrain the parameter space. 
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FIG. 2: The solid line shows the actual dependence of particle 
velocity as a function of a = a/a/ for (rn^ = 2.5keV,Xf = 
2, ^ = 0.3) while the dashed line shows the approximation we 
have used for our numerical code. 



IV. BOUND FROM FREE-STREAMING 

Generally, hidden-sector dark matter does not need to 
freeze-out when nonrelativistic (cold) in order to match 
observations [TU]. The dark- matter particles may have 
significant thermal momentum at freeze-out (it can be 
warm or hot from the hidden-sector perspective), and 
this can lead to a significant free-streaming length that 
can suppress the matter power spectrum on scales as 
large as those corresponding to galactic length scales. 
While the freeze-out process is thermal in the hidden sec- 
tor, since the hidden sector is at a different temperature 
than the visible sector, it mimics a non-thermal freeze- 
out process. While the detailed particle physics model 
can be rather different, it is quite similar to the familiar 
case of sterile neutrino warm dark matter |17H25| with 
the important difference that it did in fact thermalize in 
the hidden sector at a temperature lower than the visible 
universe. We will see that this difference will make the 
free-streaming bound somewhat less stringent compared 
to the a standard warm dark matter case with ^ = 1. 

It was pointed out in Ref. that in the case of 
warm dark matter, a common approximation for the 
free-streaming scale used for standard model neutrinos, 
Xfs{t) — 2n/kfs{t), does not correspond to the actual 
suppression scale in the matter power spectra. Instead, 
the largest scale affected due to free-streaming is essen- 
tially the particle horizon of dark matter particles. The 
comoving free-streaming horizon of a dark matter parti- 
cle with typical velocity (w) is given by 



kFSH 



-dt 



(14) 



calculation shows the momentum average of the particle 
velocity, (u) = {p/ \/m? + p^), is given by 



(v) = («(a)) 



/(9,r'') 



(15) 



9 ^ 
and 
at 
a 



Pf = 

rph 

freeze-out. 

a/af, 



dp is the momentum at freeze-out, 
(Tf is the hidden-sector temperature 
Here p is the physical momentum, 
af is the scale factor at freeze-out, 

is the Fermi- 



and /(g,rj') = [1 + exp(v/^^?T^/Tj«)]-i 
Dirac distribution at the hidden-sector freeze-out tem- 
perature. At early times a — > and, provided the freeze- 
out temperature is high enough, (u) — > 1 as expected. 
But as a increases and the momentum of each particle 
decreases we find (v) cx a~^. 



While Eq. ( 15 ) quantifies exactly how to calculate the 



thermal average, a simple approximation can be used to 
calculating A^'^^ to sufficient accuracy. We split the free- 
streaming epoch into two regimes by defining a transition 
scale factor a„,- when dark matter particles become non- 
relativistic after freeze-out at a/, and write X^^^ as a 
sum of two integrals 



kFSH 



1 



da - 



. da 



(16) 



Within our approximation, prior to a„r we take {v) 



while for a 

max[a/, a„rl 



we take (v) 



cx a 



and a\ 



Here, anr corresponds to a scale fac- 
tor when dark matter particles becomes nonrelativistic 
(T,^^ ~ m^/3.15) [To]. This corresponds to a visible 
sector temperature Tnr — {m^/ {3.15£,)){gj / g^^Y^^ . For 
particles that freeze out when already nonrelativistic, so 
that a/ > anr, only the second term contributes. In Fig- 
ure [2j we plot the numerical solution of free-streaming 
velocity as a function of scale factor and also show 
the approximation we have adopted. We find the free- 
streaming length computed using this approximation ver- 
sus Eq. ( 15 1 agree at the few percent level — accurate 



enough for this work. 

For relativistic decoupling it is straightforward to show 
that the free-streaming length is only a function of mass, 
but for the semirelativistic and nonrelativistic cases it 
depends on both mass and ^. In all cases the co- 
moving free-streaming length must satisfy the bound 
from the measurements of the linear power spectra on 
small scales which in turn constrains the dark-matter 
parameter space (^, cr, m^). We use here the bound 
^FSH < 230 kpc [inillS]. 



where tf and teq corresponds to time of freeze-out and 
matter-radiation equality respectively and a is the scale 
factor. The scale factor dependence of the velocity 
changes with Xf and the mass of dark matter. A simple 



V. TREMAINE-GUNN BOUND 

A robust and model independent lower bound on the 
mass of dark matter particles is obtained by bounding 
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FIG. 3: Representative constraints in (cr/ao,^) plane for = 1900 keV (left panel), 3 keV (middle panel), and 12.5 keV 
(right panel). Here, ctq ~ GeV~^. The horizontal dashed lines are contours of free-streaming length while the bold one 

corresponds to X^^^ ~ 230 kpc. The solid lines are contours of minimum Tremaine-Gunn mass, m^^„(^,(T), and the bold 
solid line (right panel) corresponds to m^'f„(^, a) — m^. We see that the free-streaming bound is more stringent for low dark 
matter mass, while the Tremaine-Gunn bound restricts the parameter space for higher (e.g., right panel). The colored 
region corresponds to the allowed values of dark matter density (0.1 < QoAih^ < 0.114) and that also obey the Free-streaming 
and Tremaine-Gunn bounds. The corresponding values of hidden-sector Xj = rn/Tf are shown in the sidebar to illustrate the 
nature of decoupling (i.e., whether relativistic, nonrelativistic or semirelativistic). We see that for low mass is smaller and 
decoupling tends to be relativistic, while for higher mass x'j- is larger with values corresponding to either semirelativistic (warm) 
or nonrelativistic freeze out. 



the phase-space density evolution of small galaxies like 
the dwarf spheroidal satellites (dSphs) of the Milky Way. 
Originally suggested by Tremaine and Gunn in Ref. [27] , 
it has been developed further by several authors in, for 
example, Refs. [25ti31j . If the dark matter is fermion it 
gets a stringent bound, independent of cosmological evo- 
lution, as Pauli blocking enforces a densest packing of the 
dark matter phase space distribution. However, a differ- 
ent bound relevant for both bosons and fermions applies 
to any bath of particles once in thermal equilibrium with 
interactions that freeze out. Since the microscopic phase 
space density (PSD) is exactly conserved for coUisionless 
and dissapationless particles by Liouville's theorem the 
coarse-grained (averaged) PSD can not be an increas- 
ing function of time. This means that the maximum 
coarse-grained PSD in a galaxy today, with core radius 
Tc and velocity dispersion , must not exceed the maxi- 
mum value of the microscopic PSD (the Fermi-Dirac for 
a fermion) over cosmic history f^axi^) — fmax- This 
inequality translates to a bound on the mass of 

( 9h^ 1 

(17) 

where f{q,T^) is defined in the previous section. Here 
one assumes that the collapse of dark matter from the 
initial state to bound halos is dissipationless and coUi- 
sionless. The maximum of of the distribution /(g,Tj) 
is a function of Xy, and through this the Tremain-Gunn 
mass rrijnin also depends on ^, and a. There has been 
many observational studies employing the Tremaine- 



Gunn lower bound on dark matter mass. As a repre- 
sentative example we use the results of the ultra-faint 
dwarf spheroidal galaxy Leo IV [32l [33] to restrict our 
parameter space. 



VI. RESULTS 

We now discuss our results for the allowed parameter 
space in the rather general hidden dark matter models 
discussed in this work. We performed a survey of the pa- 
rameter space in the (^, a) plane for fixed and in the 
(^,m^) plane for fixed cr — solving the freeze-out equa- 
tions numerically. We then require that hidden-sector 
dark matter: i) accounts for the observed dark-matter 
density; ii) has a free-streaming length consistent with 
structure formation; and iii) satisfies the Tremaine-Gunn 
PSD bound. Our results extends the work of Ref. [8] for 
several reasons. Firstly, we allow the size of the dark 
matter scattering cross section to vary. Secondly, freeze- 
out and decoupling does not necessarily occur when the 
dark matter particles are nonrelativistic (xy > 3). Fi- 
nally, since these constraints are relevant for relativistic 
and semirelativistic models, we include constraints from 
the free-streaming length and PSD. A special case of our 
study with tr = aEw ~ 90cro and a;^ > 3 corresponds to 
WIMPless dark matter scenario studied in Ref. [5], and 
we have checked that our calculation reproduces their 
results in that limit. 

In Fig. 3 wc show the allowed parameter space in the 
(^, <j) plane for different choices of dark-matter mass. We 
find that for dark matter masses < 1.5 keV there are 
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FIG. 4: Representative constraints in the (m^,^) plane for ct/cto = 10"^, 10~^ and 10"^. Where tJo = 10~® GeV"^. The 
colored region corresponds to the allowed values of dark matter density (0.1 < floAih^ < 0.114) that also obey free-streaming 
and Tremaine-Gunn bounds. The corresponding values of the hidden-sector x'f = m/T^ has been shown to illustrate the 
nature of decoupling (i.e., whether relativistic, nonrelativistic or semirelativistic). We see for smaller mass freeze-out tends to 
be relativistic and for larger mass it is semirelativistic or nonrelativistic. 



no values of ^ or a that satisfy all the three constraints. 
This puts a lower bound > 1.5 keV on the mass of 
hidden-sector dark matter that can account for all dark- 
matter in the Universe in a model-independent way. This 
bound is weaker than the limit on conventional = 1) 
wark darm matter (WDM) of m > 2 keV described in 
Refs. [T71 [SH [3S] because the hidden sector is at a lower 
temperature (^ < 1). 

It is clear from Fig. 3 that at very low mass the free- 
streaming bound is more important than the Tremaine- 
Gunn bound, and so the lower limit on mass > 
1.5 keV is determined mainly by the free-streaming 
bound. At higher masses the Tremaine-Gunn bound be- 
comes relevant as free-streaming becomes more and more 
suppressed (see the right panel of Fig. 3) and the al- 
lowed region of parameter space increases. Eventually, 
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FIG. 5: The allowed region in the (m^,^) plane for 
a — 1.5 X 10~® GeV"'^. As this cross section is of the same 
order of a canonical cold WIMP, the allowed region extends 
all the way up to ^ = 1 as expected for ~ 100 GeV. 



for > 25 keV, only the constraint from dark matter 
relic density remains relevant. The bands shown repre- 
sents the width arising from the range of allowed dark 
matter relic density 0.1 < fldh'^ < 0.114. The legend 
color bar represents different values of rrij^/T^ at the 
time of freeze-out and shows the transition from rela- 
tivistic, to semirelativistic, to nonrelativistic freeze-out 
as increases. We see that freeze-out can be easily 
semirelativistic 1 ^ a;/ ^ 3 or relativistic 2:^ < 1 and thus 
showing hidden-sector dark matter that freezes out when 
warm or hot is consistent with all cosmological bounds 
and may be considered when building particle-physics 
models of dark matter. 

In Fig. [4] we show the allowed parameter space in 
(^, m^) plane for three choices of the scattering cross sec- 
tion a. We see that, intuitively, for lower cross sections 
and interaction rates the freeze-out process tends to be 
relativistic, but as the cross section increases freeze-out 
can occur later and be semirelativistic or nonrelativis- 
tic. In Fig. [5] we show that we recover the standard cold 
WIMP scenario for higher mass (~100 GeV) and typical 
wimp scale cross section. As we move to higher mass, ^ 
need to be below unity to satisfy the dark matter abun- 
dance. 



VII. NONGRAVITATIONAL DETECTION OF 
HIDDEN DARK MATTER? 

If the dark matter is in a hidden sector then it has at 
most feeble interactions with the standard model so that 
it never equilibrated with the standard model plasma (or 
at least fell out of equilibrium very early on). Its direct 
detection is extremely unlikely and one would expect its 
detection only through astrophysical signatures. How- 
ever, a more precise statement is that any possible inter- 
action rate with the standard model must be less than 
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the Hubble rate TjlT) < H{T) throughout cosmic his- 
tory. In spite of this constraint, it is not impossible to get 
direct or indirect detection of hidden-sector dark matter. 
This occurs because of extra freedom in the hidden sec- 
tor freeze-out process. In contrast with typical WIMP 
models, a wide range of dark matter mass can produce 
the observed thermal relic density, and importantly the 
dark matter annihilation cross section a determines only 
the thermal relic density while the direct or indirect de- 
tection rates are determined by aj. In the context of the 
WIMPless scenario Ref. [3S| has suggested that scalar 
hidden dark matter may have non-gauge interaction with 
the standard model through connector particles, with 
interaction of the form Cx = XjXYlJl + XjXYrJh, 
where X is a scalar hidden dark matter particle, 1^,^ 
are chiral fermions and fh.R are standard model fermions. 
With such interactions hidden-sector dark matter can al- 
low interesting direct and indirect detection phenomenol- 
ogy 36 . For a fermionic hidden dark matter detection 
might be more difhcult as low-energy interactions involve 
higher dimensional operators such as the dimension 6 
four-fermion interaction (xx)(//)- This type of inter- 
action scales like cr/ ~ A^/s above an interaction scale 
TO/ and as aj ~ (A|/m|)s far below it, where s ~ ^'cm- 
Demanding that the hidden sector remain out of ther- 
mal contact with the visible sector likely requires that 
this effective interaction is extra-weak [10, (weaker than 
the Fermi interaction strength Gp), but detailed pre- 
dictions for nongravitational signatures of hidden-sector 
dark matter are best done on a case by case basis. 



details of hidden-sector particle physics and relies on 
gravitationally-mediated signatures of dark matter from 
cosmology. 
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Appendix A: Surface of Allowed DM Abundance 

We have determined a fitting formula which might 
be useful for characterizing thermal hidden dark mat- 
ter models of the kind we study here. We provide a 
polynomial fitting function for hidden dark matter scat- 
tering cross-section cr^ = fi^x^O which gives the ob- 
served dark matter relic density. This formula is valid 
for a general thermal warm or cold hidden dark matter 
model with a dark matter mass range keV to 100 

GeV, range for the hidden-to-visible sector temperature 
ratio from ^ ~ to 1, and for cross-sections ranging from 
10-3 CTo to lOo-Q where CTq = 10-^GeV-^ This function 
can be used to estimate for a given dark matter mass 
TO^ and hidden-to-visible temperature ratio ^. 



VIII. CONCLUSION 



Dark matter might belong to a hidden sector with 
new particle physics totally unknown to us. Even if 
this is the case, we can put constraints on its particle 
properties from its gravitational effects on the Universe. 
The constraints we discuss here are valid for a wide 
class of hidden-sector models of dark matter where the 
dark matter is a thermal relic particle that was once 
in equilibrium with a hidden-sector plasma. We have 
allowed the dark-matter annihilation cross section to be 
a free parameter and the hidden sector temperature to 
be different from the visible sector. Due to these extra 
freedoms we have shown, in contrast to a visible-sector 
WIMP, we must include cases where the dark matter 
is relativistic, semirelativisitc, and nonrelativistic at 
freeze-out. By solving the general freeze-out scenario 
numerically we have treated all cases in a unified 
way, and have found the region of hidden dark matter 
parameter space that can account for the observed 
density of dark matter while remaining consistent with 
current constraints on the free-streaming length and 
phase-space density of dark matter. A thermal relic in 
a hidden sector with dark matter mass < 1.5 keV 
is incompatible with current constraints. This lower 
bound on the mass of hidden-sector is insensitive to the 




40 60 
(GeV) 



I 



FIG. 6: Density plot of Eq. (All showing (ct^/ctq) values in 



the (^, m^) plane that yield the allowed dark-matter relic den- 
sity. The legend maps the value of (a^/ao) to an intensity or 
color. The grey area indicates a region where the cross section 
is less than 10"^ ctq 
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i.j=0 

where the mass is written in units of GeV and ctq — 
10~^GeV~^. We have checked that the above fitting 
function is accurate to within 4 or 5 percent by compar- 
ing it with the actual solution of the Boltzmann equation. 



o 


I = U 




t = Z 


1 = 6 


j = 


-5.2345 


-0.0016 


0.0001 


1.028 X 10"" 


j = i 


23.3984 


0.0053 


-0.0004 


-4.0929 X 10"" 


i = 2 


-29.4988 


0.0029 


0.0005 


5.4513 X 10"" 


i = 3 


15.2099 


-0.0018 


-0.0002 


-2.6184 X 10"" 



TABLE L Coefficients dj of the fitting function for the hid- 
den dark matter annihilation cross section surface with al- 
lowed relic abundance. 
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